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We reexamine the possibility of employing the viscosity over entropy density ratio as a diagnostic 
tool to identify a phase transition in hadron physics to the strongly coupled quark-gluon plasma and 
other circumstances where direct measurement of the order parameter or the free energy may be 
difficult. It has been conjectured that the minimum of rj/s does indeed occur at the phase transition. 
We now make a careful assessment in a controlled theoretical framework, the Linear Sigma Model at 
large- A^, and indeed find that the minimum of r]/ s occurs near the second-order phase transition of 
the model due to the rapid variation of the order parameter (here the a vacuum expectation value) 
at a temperature slightly smaller than the critical one. 



I. INTRODUCTION 



The viscosity over entropy density ratio rj/s has recently received much attention due to the conjectured universal 
bound [l[ r\j s > for fluids describable by any quantum field theory. No experimental exception to this rule has been 
found to date, and many have been examined, such as ordinary gases and liquids like argon [2J, undopcd graphcnc 
Q, fermions near the unitarity limit 0L and more interestingly for our purposes, the hadron phases formed after 
Relativistic Heavy Ion Collisions 0, HIM Q • The existence of this bound can be inferred from dimensional arguments 
for a quantum fluid alone [h| , but the 1 /Air numerical factor has only been obtained in the context of conformal field 
theories with a gravity dual. In [ill . \v\ some possible violations to this bound are proposed. 

Important insight came from [6( as it was further conjectured that the minimum of r\J s for several ordinary fluids, 
but also for hadron matter, coincides with the phase transition. We 0, [l3[ and others [3, [H[ have gathered further 
empirical and theoretical evidence in favor of this concept. 

Its use is quite evident for the experimental program of FAIR facility (for example the Compressed Baryonic Matter 
experiment) or RHIC if run at a lower energy. The idea is simply to use the minimum of rj/s, that might be accessible 
by studying the momentum distribution of final state pions and other particles, through elliptic flow 0, [l(| for 
example, as a diagnostic to locate the phase transition and possibly the critical point in QCD. 

For small baryon chemical potential, and based on lattice QCD calculations, one expects a crossover. Therefore the 
ratio rj/s would show a soft minimum when the quark-ghion plasma is about to appear. Our theoretical estimate of 
rj/s in this case is shown in Fig.[Tl an update of that in [2|. The hadron side (mainly containing pions but also kaons, 
r\ mesons and their dynamical, clastic low-energy resonances in minor concentration) is calculated by means of the 
Boltzmann-Uehling-Uhlenbeck (BUU) transport equation and using Chiral Perturbation Theory for the low-energy 
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FIG. 1: The phase transition between a hadron phase (left curves) and a strongly coupled quark-gluon plasma (right curves) 
points out to a minimum or a discontinuity depending on the order of the transition. At jis ~0a soft minimum is expected as 
hinted in our prediction for rj/s in the hadron side. The QGP estimations are based on the calculations in @, E^i implementing 
the formulae rj/s = 5.119/ [g 4 In (2.414 g' 1 )} for N f = 3 and rj/s = 5.328/ [g 4 In (2.558 g' 1 )] for Nf = 2, with g as a function 
of temperature given in [f|. Finally, the calculation of Demir and Bass with a microscopic transport model is also given [jj. 
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interactions with the Inverse Amplitude Method for the meson amplitude unitarization. We have not included a finite 
nucleon density in our work, but this has been estimated independently in Ref. [l8j . 

For the quark-gluon plasma side we took the perturbative estimate of [H, [T3| for massless quarks with Nf = 2 and 
Nf = 3. In [2| we showed that the corrections in this side due to the finite quark masses are very small by using some 
results in [19j . 

In this contribution we are going to sharpen the statement somewhat, and show quite clearly that the minimum of 
rj/s occurs not exactly at the phase transition, but somewhat before, where the condensate has rapid variation. 

Computation of transport coefficients in a pion gas requires knowledge of the 7T7r-scattcring amplitude. Although 
in our recent work Q we have employed Chiral Perturbation Theory and Unitarized Chiral Perturbation Theory, that 
is the correct low-energy limit of Quantum Chromodynamics, the phase transition is not accessible from it, and one 
needs to resort to a separate theoretical approach for higher energies (as in Fig. [I]). This however would obscure our 
purpose of showing the correlation between the minimum of rj/s and the phase transition. 

Since some level of modeling is necessary anyway, we adopt from the start of this article the Linear Sigma Model 
(LctM) in the limit of large-A, that although less precise to provide us the scattering amplitudes necessary in the 
computation of transport coefficients (see Appendix [B| , it allows us to control the effective potential and phase 
transition. Thus, the different phases of the model can be obtained -as desired- from the same partition function. 

An initial study in this direction is [bH ]. However, we feel that we can improve and better motivate the char- 
acterization of the phase diagram and the location of the critical temperature in the system without resorting to 
arguments of "naturalness" , through some additional mathematical effort. In this article we will therefore describe a 
more exhaustive and consistent study for both the phase diagram of the model and the calculation of rj/s. 

The rest of the paper is organized as follows. First we present minimum discussion on the LctM effective potential 
in Sec. HH but all mathematical detail is relegated to Appendix |A"1 Section ITTT1 is a brief indication of how the viscosity 
over entropy ratio is obtained in the LctM, with necessary detail in Appendices [Bl [C] and |Dj Our main point, the 
connection between the fast-changing condensate and the minimum of rj/s, is presented in Sec. IIV1 We also comment 
on the possibility of identifying the critical exponent of the viscosity near the critical end-point in the phase diagram 
in Sec. |V] and expose the behavior of rj/s in helium-4 at low temperature. Finally, Sec. I VII summarizes our main 
results and looks ahead to what interesting future work there might be. Appendices [E] and [F] provide additional 
technical aspects supporting our results, such as why the minimum of rj/s remains below T c , and the applicability of 
the transport equation. 



II. EFFECTIVE POTENTIAL AND PHASE TRANSITION 



In this Section we will obtain the phase diagram in order to characterize the phase transition temperature of 
the LctM at large-iV. The physics of the LctM at large-iV has been extensively described at zero-temperature in 



2C , |21| . |22j (for both zero and nonzero pion mass) , at finite temperature in [23j and the transport coefficients in 
2J,[25|], for instance. 
The well-known renormalizable Lagrangian density for an N + 1 multiplet <E> is 

C[^,d^} = - 9 M $ T 9^$ + 7i 2 ($ T $) - A($ T $) 2 (1) 

where ~p 2 (not to be confused with the renormalization scale ji, to be introduced in Eq. ()A13[) ) is positive (opposite 
in sign to a scalar field mass term) and A > 0. With this choice of parameters the LctM presents spontaneous 
symmetry breaking from SO(N + 1) to SO(N). The field $ acquires a vacuum expectation value (VEV) where the 
field configuration of minimum energy verifies (at tree level) 

= 2A = f * = NF2 - (2) 

Denoting 7r a (a = 1, N) to the N first components of <& and ct the N + 1 component we can choose the VEV in the 
direction of the latter. Thus we have {n a ) = but (ct(T = 0)) = f^. The pions are the N massless Nambu-Goldstone 
bosons; on the other hand, the field ct acquires a mass equal to m 2 = 8XNF 2 . Taking the limit m 2 — ► oo one can 
express ct in terms of the pions as ct = \J ] 2 — 7r a 7r a . This is the nonlinear sigma model in which one can eliminate 
explicitly the ct degree of freedom. 

Moreover, we can insert a possible 7r-mass term that explicitly breaks the SO(N + I) symmetry in the direction 
of ct. This term should be in the form £§b = cct, with c a dimensionful constant. Expanding ct in terms of the pions 
and requiring the appropriate form of the pion-mass term one can deduce that c = to 2 /^. Explicitly, 



m 2 / 2 A /l-^f=m 2 / 2 -im 2 7rV + .-. (3) 
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In the rest of the paper we will assume F is a constant and denote the pion decay constant or the a VEV at 
T = as fn (it implicitly depends on N because / w = s/NF). At finite temperature we will employ v(T) = (cr(T)) 
indistinctly. In particular, note that v 2 (Q) = NF 2 = f 2 . 

The phase diagram can be understood by computing the effective potential V e s at finite temperature as a function of 
the cr-condensate, to ascertain whether and when the minimum of the potential corresponds to a symmetry-breaking 
phase. 

This effective potential is sometimes extracted from the generating functional of 1PI, n-point Green functions, 
r(</>, G _1 ) (see [26[ for details). Here, we follow a different (but related) approach. The key formulas for understanding 
the procedure are given in Appendix [AJ 

Once computed, the effective potential is a function of the er-field, an auxiliary inverse propagator G _1 [0, X ] and 
the temperature. It reads, in terms of renormalized quantities, 



v. = \{" 2 - »n <r- M -M (i + £ 1O6 ^) + M#_ m ^. f! , 0(T , G - 1(0 , x , 

(4) 

Assuming a spatially homogeneous a condensate, the saddle-point approximation for this effective potential yields 
the value of (a) , together with a second equation of motion which is just the constraint necessary to solve for G _1 [0, x] > 

dV eS _ dV cS _ 

da ~ U ' d(G-![0,x]) ' [) 

where the analogy with the Cornwall- Jackiw-Tomboulis formalism is evident at this point (compare with Eqs. (1.1a- 
1.1b) in [Hj|). Explicitly these equations read 



Q = aG- 1 [0 lX }-mif 7T , 

= I (a 2 — 7VF 2 ) - (± - ^bg^M) + f 9l (T,G-i[0, X ]). W 



The first one minimizes the potential in terms of a. The second is the implicit constraint for G _1 [0, x] f° r an Y value 
of V e g, not necessarily at the minimum. When = the first equation in © reveals that there exist one phase with 
a = and G _1 [0, x) 7^ and another phase where G _1 [0,x] = and a ^ 0. In the latter case, the second equation 
yields 

a 2 (T)=AF 2 -A 9l (T,0), (7) 
and finally using the result in (|A20j) one arrives at 

T 2\ 1/2 



*(T) = W ( 1 - 55 ) , (8) 
with a 2 = f 2 and the critical temperature 

T 2 = 12F 2 . (9) 

One can also check that the critical exponent for the order parameter is equal to the one from mean-field theory 
13 = 1/2. 

The case with 7^ is not so easily tractable and must be solved on a computer by an iterative method, such as 
Newton's method. Convergence is easily achieved. 

Eliminating G _1 [0,x] this way, one can easily plot the effective potential as a function of a at fixed temperature 
(Fig. [2]). Identifying the position of its minimum with the help of a computer, we can also plot the value of the a 
condensate at the minimum as a function of T (see Fig. [3]). 

In Fig. [21 we obtain the shape of the effective potential as a function of the relevant parameters, temperature T, 
renormalized mass of the <r, Mr 1 and the physical mass of the pions (analogous to the coupling of an external 
field that explicitly breaks the original symmetry). For typical values of these parameters see Table |H 



1 Mr can be traded for the renormalized coupling constant A^j, through Eq. (IB5l t 
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FIG. 2: Effective potential in various regimes of the LctM, showing the various phases and the nature of the phase transitions 
between them. Top row: massless pions. Bottom row: pion mass fixed at 139 MeV. From left to right, the renormalized a 
mass and temperature are varied as indicated in the graph. Each potential has been multiplied by an arbitrary temperature- 
independent constant factor, to match their scales for visibility. See the main text for interpretation of the results. 




FIG. 3: The value of the a condensate as a function of temperature. Setting — 0, the model presents a clear second-order 
phase transition (PT) with a discontinuity in the derivative. For finite pion masses, we have a crossover. 



TABLE I: Default parameters used in the computation of rj/s in the LctM (used unless noted otherwise). 

Parameter Value 

139.57 MeV 

U 93 MeV 

TV 3 

M R0 0.5 GeV 

X R M 2 m /{%fl) ~ 3.6 

[i 2 1 GeV 2 " 

T c 2V3F = 0.186 GeV 



"Nothing depends on this scale choice when using (IB3I > and ( IBS t ■ 
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Inspection of the top row in Fig. [2] reveals a second-order phase transition from a potential with a single minimum 
at a = at high temperature, to a phase with spontaneously broken symmetry below the critical temperature, where 



An imaginary potential Im V c s ^ arises when the real part V e g becomes convex, with negative second derivative. 
This is a well-understood phenomenon [27j and is related with the probability of decay per unit volume of the unstable 
vacuum. 

In the bottom row we consider the effect of a finite pion mass. Then the potential loses the residual left-right 
symmetry (reflection respect to the a = axis), and the pion mass is acting as an external magnetic field in an 
analogous condensed matter system. The second-order phase transition becomes then a smooth crossover. 

In order to further characterize the phase transition we now turn to Fig. [3J with the dependence of a taken at the 
minimum of V e s, as a function of T. This is nothing but the order parameter of the system. The two possibilities 
in this model are now clear: a second-order phase transition and a crossover, if we introduce an explicit pion mass. 
Note that only the positive VEV is depicted in Fig. [3j 

The phase transition in this model is a chiral one, where the pion acquires a thermal mass above the critical 
temperature. This continuously matches to zero at the critical temperature, and we cannot show a first order phase 
transition for any combination of the parameters 2 . 

Combining the knowledge of V e s(cr) and (er(T)) we arc able to determine the critical temperature of the chiral 
transition and more generally, the entire phase diagram of the theory. 



In this section we quote our computations of rj/s for the LcM at large- TV. To compute the shear viscosity we 
employ two different methods. 

The first calculation uses the quantum transport equation of traditional kinetic theory. Many details for this 
Boltzmann-formulation are given in our other publications p8l . [2^ | and we omit them here for shortness. In Ap- 
pendix [C] we sketch a few steps for completeness. The Boltzmann-Uehling-Uhlcnbeck is a transport equation for the 
distribution function of the gas. This equation is linearized and solved perturbatively in an appropriate polynomial 
basis, in the so-called Chapman-Enskog formalism. The shear viscosity of the gas is then easy to obtain projecting 
the solution. 

An alternative method to calculate the shear viscosity is the Green-Kubo equation. This method, based on Linear 
Response Theory, works better for a field theory that accepts a diagrammatic expansion. Some details of the Green- 
Kubo formalism are reviewed in the Appendix [Dj Consistency of both approaches provides a check of the calculation 
through kinetic theory and field theory methods. Results should be equivalent as has been shown in [3C| . 

A further check is provided in Appendix [El where we revise the validity of the Boltzmann approximation to out 
of equilibrium computations. In essence, this "molecular chaos" hypothesis needs successive collisions for the same 
particle to be uncorrelated, requiring small densities as compared to the typical scattering cross section, ncr 3 / 2 << 1. 
This we establish in that Appendix for the range of validity of our study. 

The other ingredient needed for rj/s is the entropy density, that can also be calculated in two different ways. 
Frequently quoted is the free entropy of a degenerate Bose gas. For vanishing chemical potential one can obtain a 
simple formula for the entropy density from the Bose-Einstein partition function 



(a) jk 0. 



III. COMPUTATION OF THE KSS NUMBER: v /s 




(10) 



(11) 



For massless pions, this simply reduces to 




(12) 



2 One can artificially force it by introducing an ad-hoc discontinuity in m„- in the computer code, and we have examined this case to 
check that rj/s inherits the same discontinuity. 
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FIG. 4: Comparison of the viscosity computed with the Boltzmann-Uehling-Uhlenbeck equation (left panel) and the Green- 
Kubo formalism (right panel). The result is nicely consistent for a broad range of temperatures and renormalized masses. 



An a priori better and more accurate alternative is the entropy density derived from the effective potential in 
Eq. (H]). One should then regard V e g as a Hclmhotz free energy obtained from the generating functional (as a 
partition function). 

, m _ _ dV cfi _ N dgp^G-^x]) 

[ j ~ dT ~ 2 dT ■ [U) 

The limit of Eq. (p~2|) for noninteracting massless pions is recovered upon employing Eq. (|A19[) . We again resort to 
a computer for the, more complicated general case; our numerical evaluation however shows that the entropy remains 
close to the free gas, since at moderate temperatures the density is small. 

Combining the results for rj and s we can deduce rj/s as a function of temperature and plot it in Fig. [5] Both 
formalisms are shown (left panel from Boltzmann equation, right panel from Green-Kubo formula) and one can see 
the excellent agreement between them. 

The last step in our analysis is to connect the results of rj/s and those for the critical temperature of the LerM. 



IV. MINIMUM OF rj/s AND VARIATION OF THE (a) CONDENSATE 

We now establish that the minimum of rj/s tracks the movement of the critical temperature when the model 
parameter F is varied. For this we return to the case = so as to have a second-order phase transition in the 
LaM. Recalling Eq. ([9]) the critical temperature is 

T c = 2V3F , (14) 

Therefore we need to extract the dependence of rj/s on F. This is done in the right panel of Fig. [5] For ease of 
comparison, we show in the left panel how T c depends on F. For the studied range of F the critical temperature 
(defined as the temperature at which (a) reaches a zero value), extracted from the computer code, varies exactly like 
the analytical Eq. (fl4|) . which provides a handy check. 

First, we note that rj/s presents a nonanalyticity as a function of T at T = T c where a very different qualitative 
behavior of rj/s begins. This is due to the dependence of rj/s on the thermal pion mass m^(T) and the a condensate, 
both nonanalytical at T c . Below T c , the thermal pion mass is identically zero where Goldstone's theorem protects 
the pion masslcssness from radiative corrections. In this phase the symmetry is still broken and the condensate (er) 
is nonzero and described by Equation l[8|). Above T c the symmetry is restored and (er) vanishes. The "classical" pion 
mass 771^ is zero as well. However, the quantum thermal corrections force a temperature-dependent thermal mass, 
m ff (T) that is obtained from the thermal tadpole in the pion self-energy. Thus we have 

(cr(T)) = I CT(T) T<Tc (15) 
I T > T c 



ml(T) = { I, /ma ^ t:t C (16) 
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FIG. 5: We follow the evolution with F of the minimum of rj/s and the a VEV, that yields the critical temperature at vanishing 
condensate. Both increase the same relative amount when varying the coupling constant. 




T (GeV) 

FIG. 6: Our key result: the minimum of r)/s occurs just before the critical temperature for the phase transition in the Linear 
Sigma Model. This is where the condensate varies rapidly approaching zero. The phase transition is marked by a nonanalyticity 
of the ratio. 

that is continuous at T c but nonanalytic. The behavior of (a) itself is also continuous as shown in Fig. [3] but with a 
discontinuous derivative at T c . These two quantities influence both rj and s and their nonanaliticity is inherited by 
the KSS ratio. 

By varying F, the minimum of ?y/s moves according to the F dependence of T c , providing evidence for our claim 
that T c and the minimum of rj/s are related. 

In Fig. [6] we zoom in the r//s plot near the critical temperature. As can be seen, the minimum is not reached at the 
critical temperature, but right before T c . This shows that the minimum of rj/s is controlled by the rapid variation 
of the order parameter. To understand this result one needs to keep in mind the diffusive nature of the transport in 
a gas. With increasing temperature, the gas particles carry transverse momentum between different parts of the gas 
more efficiently, and thus increase the shear viscosity. However their interactions hamper transport. This can be seen 
from Eq. (|B7[) . As (a) decreases rapidly, the pion elastic cross section increases. Since rj oc l/ow in kinetic theory, 
the viscosity must drop. Eventually growth is regained as the temperature increases. 

Although the LcrM docs not present a first order phase transition, within our treatment one can artificially force it 
by adding a jump for (a) in the computer code. In this case, the viscosity over entropy density ratio also turns out 
to have a discontinuity in the program (not shown since it is not a genuine model prediction) with a minimum of rj/s 
at exactly T c due to the jump, analogous to the first order transition in atomic Argon in our prior work (2j. 
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V. CRITICAL EXPONENT FOR rj/s 

We have argued [l3| that r\j s has a tiny critical exponent of about y ~ 0.04 at the critical end-point of QCD by 
applying the theory of dynamic critical phenomena [3lJ. This is true only if the dynamical universality class of the 
QCD critical point is that of Model H. as claimed in [33|. The critical exponent y (sometimes called <p) is defined 
near the critical temperature as 

»7 0c|i|-», (17) 

with t = (T — T c )/T c being the reduced temperature. The shear viscosity and the correlation length £ are related by 
a further power-law with exponent x v 

ri oc£ x i. (18) 

Because £ oc \t\~ u one can relate both dynamical exponents through the static one y = x^v. 

The critical exponent x v appears to be too small (near 0.06) to be detected experimentally in a collision of heavy 
ions, a rather noisy environment. For this reason we expect that experiments will bear a residual minimum (as a 
function of temperature) in rj/s, near the critical end-point. 

This does not mean that the critical behavior is unobservable. On the other hand, once measurements of the volume 
viscosity become available, its own critical exponent, near to 3 (33l |. might well be observable. 

Another not-so-small critical exponent to detect is z, defined as r a f , where r is the relaxation time and £ 
the correlation length. Near the critical point, £ diverges and the previous formula reflects the so-called "critical 
slowing-down" . 

It is interesting to note that the hyperscaling expression 3 in [3l], [HJ relates the critical exponents for the shear 
viscosity and heat conductivity x n + x\ = 4 — d — r). Further employing the relation z = 4 — r\ — x\ one obtains a 
relation between dynamical critical exponents: 

x v = z — d, (19) 

where d is the spatial dimension. For the QCD critical point, belonging to Model H one sees that z > 3 and x v 
is nearly zero. An estimate of x^ might be obtained with the assistance of Eq. (|19|) . but needing a very precise 
measurement of z through the long-time behavior of the relaxation process. 

It is highly possible that a few modes alone are responsible for this theorized divergent behavior of the viscosity 
at the critical end-point, while other system modes yield simply a broad minimum. This separation is sometimes 
referred to [34[ as "regular viscosity" and "singular viscosity". In any case, from the LoM side, mean-field theory 
(sort of what one recovers in large N) loses all information about the critical exponent (remember for example the 
3D Ising model where the critical exponent for the heat capacity is a ~ 0.11 and mean-field theory yields zero). 

In spite of the presumed critical behavior of shear viscosity in systems belonging to Model H (such as, presumably, 
QCD, but not the LoM here treated), there exist many systems where r\ does not diverge at the critical point. The 
divergence of rj/s is therefore not truly universal at the critical point, it depends on which universality class the 
system belongs to. The presence or absence of this and other critical exponents is important to help us classify the 
universality class and dynamical model of the QCD phase transition. 

A. The case of Helium-4 

As an example, we will bring forth the supcrfluid phase transition of Helium-4 at the A-point. Its dynamic 
universality class is that of Model F, in which the order parameter is taken to be the macroscopic wave function, ip. 
The quantity n = |?/>| 2 tells about the concentration of quasiparticles belonging to the condensate. In this universality 
class [3lj 77 is not singular at T c . Without passing by the critical end-point, the behavior of r\j s as a function of 
temperature is shown in Fig. [7J The data are experimental measurements for gaseous Helium, normal liquid He-I and 
supcrfluid He-II. We first point out that the KSS bound is always preserved, even at T\. Below T\ the superfluid is 
formed by the condensation of Bogoliubov quasi-particlcs. The nonzero value of 77 and s in He-II can be understood 
with Landau's "two-fluid" model. Both coefficients approach zero when T descends towards absolute zero, but the 
behavior of ?]/s at low temperatures is to increase as observed early-on [l[. 



3 We thank Victor Martin Mayor for the observation. 



9 



12 
8 



3 



5 

CD 



^— 4 
S5- 2 



0. 



1.8 2 2.2 2.4 2.6 



Helium-4 @ P=0.01 MPa 
KSS Bounds 1/(4*) 



2 3 4 

Temperature (K) 



FIG. 7: r?/s for the three phases of helium-4 at low pressure, namely, He-II, He-I and gaseous helium. 



More interesting is that two apparent local minima occur at both critical temperatures (supcrfiuid T\ and vapor- 
ization T c ). A double minimum would entail that rj/s could have a maximum between T\ — 2.177 K and T c = 2.489 
K, in the He-I phase. (This phase is not a conventional liquid, the viscosity itself is anomalous and it softly decreases 
when temperature decreases). Due to the low resolution of the experimental data, clear minima are not seen, it is 
interesting to conjecture whether this behavior could be truly universal, i.e., independent of the dynamical universality 
class, with more generality than our theoretical considerations based on the a model. This is left for future work. 



VI. SUMMARY AND OUTLOOK 



In this paper we have been able to obtain the entire phase diagram as function of (T, Mr, m ff ) of the LcM in the 
large-iV limit by means of the effective potential of the theory. We have located the critical temperature where the 
chiral phase transition takes place. The temperature dependence of the order parameter, i.e., the condensate (a) has 
also been showed. 

We have calculated the KSS number, or tj/s, for this model as a function of temperature using both quantum 
kinetic theory and the Green-Kubo formalism. A very good agreement between them has been obtained for every 
value of Mr, and other parameters of the model. 

Combining both results we have shown that a strong dependence indeed exists between the minimum of rj/s and 
the phase transition of the LoM at large-A. This minimum is not achieved exactly at T c where the phase transition 
occurs, but a little before, due to the falling of (a) towards zero. 

We also would like to discuss the generality and reach of the result. The minimum of rj/s near the phase transition 
does not depend on any parameter exclusive of the LcrM. It just depends on the falling behavior of (a) from a nonzero 
value at T = to zero at T c . The VEV of the a field can be understood as the order parameter of the chiral phase 
transition, and every other phase transitions has an equivalent order parameter (spontancus magnetization in spin 
systems, macroscopic wave function in superfluids and superconductors...). If this order parameter influences the 
transport coefficient, one should expect our conclusion to be generalized in a straightforward way. 

We have scanned the recent literature to see if other fluid systems that have been considered in connection with 
the KSS conjecture have a minimum value of r\j s that occurs at a temperature discernibly smaller than the phase 
transition. This behavior is indeed present already in but it seems to have passed unnoticed. In Fig. [5] we replot 
the data given in Q but add the location of the critical temperature claimed by the same experimental group that 
provided the original data. The minimum of r\j ' s for this cold Fermi gas near the unitary limit occurs slightly before 
the critical temperature. We look forward to other examples where enough precision can be achieved to separate the 
phase transition and the minimum of the viscosity. 

Meanwhile, one could ask what the situation for Yang-Mills theory is (as relevant to the Quark-Gluon-Plasma) . 
The perturbative computations by other authors that we have quoted [17| are not of much use near T c . Using a 
nonperturbative computation, Antonov has shown [35| that a possible minimum should occur in the range 2T C — 3T C , 
definitely above the phase transition. This could fit our observation that a significant change of the order parameter 
can mark the minimum of rj/s, as it has been pointed out by Hidaka and Pisarski [36| that the Polyakov loop acts as 
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FIG. 8: The plot from 0] shows that the minimum of rj/s for a trapped cold atomic gas in the unitary limit appears slightly 
before the phase transition temperature T c = 0.29(2)Tf. A simple quadratic fit of the data shows a minimum just below 
T = 0.2T F with a yf/d.o.f = 0.071 whereas if we constrain the minimum at T c /Tf, the quadratic fit has an almost six times 
larger x /d.o.f = 0.39. 

an order parameter and influences rj/s. The Polyakov loop only reaches a value near 1 at about 2T C . 
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APPENDIX A: THE EFFECTIVE POTENTIAL FOR THE LINEAR SIGMA MODEL AT LARGE— /V 

Taking the original Lagrangian in ([1} and performing a Wick rotation (t — > — it), the Euclidean or imaginary time 
Lagrangian is (in the following we denote <£> 2 = $ T $) 



1 



where now the product in the //-index is to be understood as Euclidean. The partition function reads 



Z = f[d$] exp(-5 B [*]) = Jldir a }[da}cxp(-S E [ir a ,<j}). 



(Al) 



(A2) 



We wish to obtain an effective potential for the a field that allows us to decide under what circumstances one is in 
a condensed phase. For this we need to perform the integration over the dynamical 7r a variables. Since the action is 
quartic in them, we introduce an auxiliary static field Xj with no kinetic term, that splits the quartic pion vertex. To 
assist large- N counting, it is convenient to include a factor of N^ 1 / 2 in its definition, 



(A3) 



The quartic potential is then hidden in the identity 



exp (-X J d 4 x ($ 2 ) 2 j = J [dx] cxp i J d 4 x (n x 2 - x$ 2 V8Aiv) 



(A4) 
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X 



FIG. 9: Trilinear coupling between <E> and the auxiliary field x that replaces the original quartic interaction $ 



that features a trilinear coupling instead and is therefore quadratic in the pion fields. This coupling is shown in Fig.[SJ 
This method has also been used in the LcrM at T = in [2fJ . 
Then, the complete partition function becomes 



Z = J {dTr a }{d<j}{d X }exp(-S cS {K a ,v,x}) 
with the full (still Euclidean) effective action 



S eS [7r a ,a,x] = J d 
Integrating by parts the 7r a kinetic term we obtain 



2^ a ) 2 + 2^) 2 - -PV - 2^X 2 + -X^VSXN + -x^VSXN - mlUa 



(A5) 



(A6) 



S cS [TT a ,<j, x ] = / d 



-7T a (-a E - 2-f + X V8XN)TT a + -(d^a) 2 -JI 2 a 2 - -N X 2 + -x^VsXN-mlUa 

To make contact with the Cornwall- Jackiw-Tomboulis formalism in [26| , it is useful to further trade x by the amputated 
two-point function for the er, that can then be used as the auxiliary variable. 



G- 1 [q, X ] = q 2 -2jI 2 +xV8XN. 
We may then perform the path integration over the 7r a fields to obtain the effective action for the a alone, 



(A7) 



[d-K a \ exp 
= I d x exp 



1 



N 



d A x ir a (-a E + G- 1 [0,x})Tr a 



J d A x [det(-D B + G- 1 [0,x])]' 



■N/2 



__ t,. log(-n B + G- 1 [0,x]) 



= Id x exp 



N f d q i n-M i 



(A8) 



Therefore the effective action to be employed in the generating functional Z = J[da][dx] exp (— r c ff[<r, x\) i s i m_ 

27rn 





troducing the finite-temperature Imaginary Time Formalism with Matsubara frequencies ui n = ((3 = 1/T) , 



f d 4 q->TJ2 n f d 3 q: 



F eS [a, X } = Jd D xU(d fi a) 2 -^x 2 + l<r 2 G- 1 [0, X ]-rnlUa+^T £ J ^logtT^x] 
From Eq. (|A9[) we read-off the effective potential. Rewriting the —Nx 2 /2 term in terms of G _1 [0,x] 



(A9) 



N X 2 _ (G- 1 [0, X ] + 2JI 2 ) 2 _ (G-MO.xl) 2 , XN 2 F , , NF^G-^x 



16A 



16A 



we obtain V e g (a, x) , 



(A10) 



V cS (a,G- 1 [q,x]) = l(<? 2 -NF 2 )G- 1 lO,x]~^(G- 1 lO,x]) 2 -\N 2 F*- m lUa+^T j[] f -j^L logG"^ x], 

n——oo ^ ' 

(All) 
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As is well known, the thermal loop integrals such as in Eq. (jAllI) contain no new divergences aside from those in 
the vacuum integral (the term n = 0). In this case 

/ dq log (q 2 + M 2 ) = T ^2 / d 3 q log ( W » + q 2 + M 2 ) = —I M - —-- - g {T; A/ 2 ), (A12) 

where the zero-frequency term in the Matsubara sum, the vacuum loop integral, carries the divergence, that one can 
extract in dimensional regularization with N e = - + log47r — 7, 



'A/ 



M 2 ( ii 



(4tt) 2 V b M 2 



/V £ + l + log^ (A13) 



and go is the thermal part of the loop integral. These p-functions are those previously defined in [37]]. Another useful 
integral (needed in the saddle-point equations) is the following: 

l d « wtw* = T E / ^ z^twtip = lM + * (T; M2) - (A14) 

We need one more thermal p-integral since Eq. ([6]) is solved by iteration using Newton's method, which requires a 
derivative respect to G -1 [0, x]- The series of ^-integrals satisfy the recursion relation 

9n+1 =-jMf (A15) 
In terms of the variable y = M/T, the ones we need are 

rpi roc 1 

50 (T;M 2 ) = d^-^-L- (A16) 

t 2 r°° , /™2 _ ,.2 
^«*)-5?/, * J Wf- (A17) 
If 00 . 1 1 



52 (T;M 2 ) = / _1_ . (A18) 

In the limit M <C T, or y — > 0, one recovers the standard results for a massless Bose gas 

rpA fOO 3 rp4 rpA A ^2^4 

*>< T ' °) = 3^ / i^Tf = 3^ 4 >« 4 > - 3^ 6 90 — ^ 



T 2 Z" 00 x T 2 T 2 7T 2 T 2 

91 (T,0) = ±^ ^ = ' m C(2) = ^ = L- (A20) 
27r z J — 1 27r z 2-k z 6 12 



As the Bose-Einstein factors cut-off high momenta, only the vacuum parts contain divergences, and these can be 
reabsorbed in the tr-mass and the coupling constant that appear in the classical Lagrangian density. In the next 
Appendix we describe the renormalization procedure. 

APPENDIX B:w-tt AMPLITUDE IN THE LctM AT LARGE-TV 

In a nutshell, transport in a gas is a diffusive effect, so that in kinetic theory transport coefficients are proportional 
to the mean free path, or inversely proportional to the total collision cross-section. Thus we need the scattering 
amplitude in the LcrM in the large-TV limit. 

This 7r — 7r scattering amplitude at tree-level is simply the one shown in Fig. 1101 



Combining both diagrams, 

Ms) = 4jtV. (Bl) 
u 1 TP 
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M") = 






FIG. 10: Tree level vr - vr amplitude at <D{\/N). 



where v is the VEV of the a field, that depends on the temperature v(T) = (<t(T)) and M 2 



$>\v is the 



spontaneously generated mass of the longitudinal mode h. This exchanged boson h is the physical quantum over the 
vacuum, a = v + h. In (|Blj) we explicitly see that Aq is 0(l/N). In the large-iV limit (with NF 2 fixed) the s-channcl 
iteration of the tree level diagram is also of order 1 /N (recall that every pion loop carries a factor N) , and must be 
resummed. The situation is depicted in Fig. 111! 



A(s] 





- 




7T 



7T 



FIG. 11: Resummed amplitude for n — n scattering at 0(1/N). 



The one-loop integral is I(s) = [7V e + 2 + Log (— [J 2 /s)~\ , where the Log is to be understood as complex and the 
integration has been regulated with the procedure of dimensional regularization. In this scheme one introduces a scale 
fi in the definition of the renormalized mass, Mn(fi) in order to absorb the infinity in I(s). 
Defining Mr as the value of the renormalized mass at some predetermined energy //q 



Mr 



Mr(ji ), 



we can express the mass at an arbitrary renormalization scale fi 

M 2 R {H) = 



Ml 



1 - 



Mi 



log 



(B2) 



(B3) 



32n 2 F 2 

Not only the cr-mass is renormalized but also the coupling constant, A. The dressed coupling constant is defined as 



1 



1 

A 



N 
4^ 



2 + log 47r - 7 + 



Both renormalized parameters are related (like their bare counterparts) through 

M R (fi) - 8X R (fi)v 2 . 

The resummed amplitude in Fig. [11] becomes 

A (s) s 1 



A(s,t,u) = A(s) 



1 



NI(s 



■Ms) 



v 2 ^ 



sN 



Log 



(B4) 



(B5) 



(B6) 



The Mandelstam variable is positive s > 0, so that choosing the branch cut of the logarithm along the positive s 
axis, in order that Log (— s) = logs + in, we arrive at 



14 



The squared scattering amplitude is simply 



\A(s)\ 2 = -jG(s), 



(B8) 



in terms of the loop function 



G(s) 



(B9) 



1 ~ Tif~(a) + 32 S X^ l0 S (p) +(3!^) 

If we consider the partial amplitudes projected on the isospin channels we see that only To dominates at order 1/N: 

(BIO) 
(BH) 
(B12) 



T = NA(s) + A(t) + A{u) = NA(s) — 



Ti = A{t) - A{u) = O 
T 2 = A(t)+A{u) = 



It is also useful to construct the isospin-spin partial waves. They are dchncd as 



t u (s) = / d( cos 9)Tj ( s, cos 9) Pj( cos 9), 

64tt /_! 



(B13) 



where Pj arc the Legendre polynomials. For moderate energies (|p| < 1 GeV) we only keep the lowest partial wave 
in each isospin channel, namely I J — 00, 11, 20. In the limit TV — > 00 their TV-dependence is 



*oo(s) 
t n (s) 



NA{s) 



O 

i 20 (s) = O 



32tt 
1 

N 
1 

N 



■oil 

1 N 



(B14) 
(B15) 
(B16) 



so too(s) dominates the amplitude (note in true two-flavor QCD at N = 3, in(s) is also very strong due to the p 
resonance) . 

For the case in which the pion mass is to be considered the amplitude A(s) must be modified. The extended 
discussion of this case can be found in [2ll ] . The partial wave is expressed as 



N 
6Zn 



(B17) 



where the amplitudes are respectively 



Aq(s) — —- — — -, 

v 1 ~ TT^JT) ~ 32^2 1 (s; H) 



,2 1 



A m {s) 



2s sN 1 



2 ' 



1-^log 



1 



4m'i 



1 



Ami , / ttvL 
1 - - log f 



(B18) 



(B19) 



When m^ = 0, the amplitude A m (s) vanishes and T(s;/i) — > log(— H 2 / s) and the (IJ) — (00) partial amplitude 
reduces to 
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l*oo| 2 (*) = T^T-^-r ~ = 7^-^ G ( s )> ( B2 °) 



(32^)2 ^ _ ^ + ^ /^^ 2 ( 32 ^ 2 ) 



, 32x1)^ , 



with G(s) defined in (|B9j) . The same result could be obtained from | too | 2 in (jB14j) . 

From Eq. (|B 14|) . the large-TV scattering amplitude |T| 2 is given by the s-channel, \T\ 2 ~ |A(s)| 2 . The averaged 
cross section in the CM frame follows 

^ _ i m 2 rR21 ^ 

dfi~ 2 647rV ( ' 

or 

= (B22) 

is the total cross section averaged over the final states. 

On the other hand, the total cross section can be defined in terms of the isospin-spin projected scattering amplitudes 

as 

3,s 

and the integrated total cross section, when only (IJ) = (00) is considered, is 

. 32tt s 2 N 2 sN 2 
CTtot(a) o /on 2^2 G 5 = a7T^ G s B24 

Note that the Kubo formalism calls for the total cross section summed over isospins, while in Boltzmann's collision 
term the initial flavors are averaged, so there is a flavor factor between them 

^ = — . (B25) 
a 3 

To compare with experimental ir — ir data at N = 3, the scalar phase-shift Sqq corresponding to Eq. (|B20j) is plotted 
in Fig. [T3J and reasonable agreement is seen. As discussed in |2l[ data are well reproduced by taking Mr — ► oo, 
the only remaining parameters being F and [i. This renormalization scale is conventionally taken to be the p mass, 
fj, ~ 775 MeV, where one naturally expects (because we have neglected the contribution of the (I J) = (11) channel) 
that the LcrM ceases to be valid for N = 3. 

An interesting point is to examine the breaking of scale invariance even at the critical temperature. This can 
be followed from Eq. (|B20[) and (|B23[) . All one needs to do is note that limr^T c V (T) — as the condensate 
vanishes, remembering also the relation of Mr and v through Eq. (|B5j) . Then the cross section fails to scale as 
a oc 1/s (conformality) due to the renormalization scale that brings in the usual logarithm of s/fi 2 . This is simply a 
consequence of scale-invariance being anomalous. 



APPENDIX C: COMPUTATION OF THE SHEAR VISCOSITY IN KINETIC THEORY 



One approach to calculate the transport coefficients is based on kinetic theory by solving a transport equation 

df(x,v,t) 



dt 



C[/(x,v,t)] 



(CI) 



that is an integro-differential equation for the one-particle distribution function /(x, v, t). This distribution function 
depends on time because we deal with a nonequilibrium gas. The explicit equation for our boson gas is called the 
Boltzmann-Uchling-Uhlcnbcck (BUU) equation and reads 



1 + 



N 



N 



1 



(C2) 
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FIG. 12: The scalar phase shift in the large-mass, large- TV limit of the LaM compared with the experimental data of [38l. l39l. l40l| . 
The massless-pion curve (dashed red on-line) corresponds to the amplitude given in Eq. (|B 18|) and the massive-pion curve is 
obtained when adding Eq. (|B18[l and Eq. (|B19j) . 



The first step is a linearization of Eq. (|C2[) . Slightly out of equilibrium, 

/^/o + <5/ = /o(l + f), (C3) 

where f is the Bose-Einstein distribution function, solution to the transport equation with vanishing collision term. 
Following the Chapman-Enskog expansion, the perturbation function \ must be proportional to the gradients of the 
hydrodynamical fields out of equilibrium. For the viscosity these are nothing but the velocity fields. Following the 
lines of [H| one can parametrize Sf as 



/ = fo 



i+ 9 M 



rp v 



where Ajj = PiPj — \5ijP 2 and Vij represents the gradient of the spatial velocity field, Vi 



1 



1 



Vij = Ti&Vj + djVi) - -d k V k 5, 



(C4) 



(C5) 



The scalar function g(p) depends only on moduli of momenta and is expanded in a convenient polynomial basis 

g{x)=Y J B s P s {x;y,z), (C6) 



s=0 



where x = p 2 /m\, y = m^/T and z is the fugacity of the boson gas. 

This polynomial functions are used as variational ansatz to systematically solve the linearized equation. 
The microscopic expression for the viscosity in terms of g(p) becomes 



(C7) 



where g{p) should be determined by solving the BUU equation. 

Taking only the first term in the expansion (|C6[) and inserting it into the viscosity we obtain 



2n Nm% „ 



(C8) 
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with the function 







and Bo = C^/A^, where 



(CIO) 



and the collision integral 

An = J daVr cl dpd Pl e^ E ~ 2m ^f foifUoi (l " ze^"^) - pV (Cll) 

p'V'(l - e - W£ '-"')-pV(l - e-^-^+p'fp'Kl - e -0(^-/O) -p k lP [(l - e -/3(£i-ri)l . (C12) 

This collision integral requires intense computational work. Once the kinematics are analyzed and conserved 
quantities employed, the collision term turns into a complicated five-dimensional integral, handled with a Monte 
Carlo integration routine. The calculation has already been described in (28j . 

APPENDIX D: SHEAR VISCOSITY IN THE GREEN-KUBO FORMALISM 

In this Appendix we briefly review the computation of the shear viscosity in the Green-Kubo formalism, based 
largely on the works of [4l|, |42[ ■ The Green-Kubo formulas for transport coefficients are based on the Linear Response 
Theory. If the response of a system to an external perturbation is considered to be linear in the perturbing field, the 
transport coefficient associated with a conserved quantity can be expressed as an expected value at equilibrium of the 
commutator of the corresponding (perturbed) Noether current. 

For the shear viscosity, 

r, = i km lim - f d 4 x e*- ([ft«(*),*«(0)]>, (Dl) 
20 w^o+ |p|->o+ uj J 1 

where tt 1 ^ = — g y T^/3. Coincidence with the transport equation approach occurs at low temperature. In this 
limit, the simplest dominant resummation [4j| gives 

,(o) i r^, ipi 6 



V = lO^T J ^WWfsiE^l + n B {E p )\, (D2) 
where E 2 = p 2 + m^, ub{E) is the Bose- Einstein occupation number 

n B {E p ) = (D3) 

and r(p) is the pion width in the thermal bath. The particle width can be understood as the inverse of the collision 
time in the gas and is related to the pion self-energy by 

m = _ nnn^ p ,p ) > (D4) 

where the retarded self-energy Hr(E p ,p) is given by the quantum fluctuations in the medium that change the 
dispersion relation of the pion. 

In the Dilute Gas Approximation ((3E P 3> 1) the thermal width is given by 



r(p) * i / q tot ( S ) Vs(S ~ 4m ' ) e"^, (D5) 

yFJ 2 J (2tt) 3 V ; 2E k E p ' V ; 



where as usual Ek = \/k 2 + m 2 , —> k in the chiral limit, s(k,p) is the Mandelstam variable and CTtot(s) is the total 
pion-pion scattering cross section for the LctM obtained in (|B24j) . 

As in the kinetic theory approach the coefficient of shear viscosity turns out to be inversely proportional to the 
cross-section. 
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FIG. 13: Left panel: Ratio of NLO to LO terms in a derivative expansion of Eq. (|B8|) . Right panel: Ratio of full to LO squared 
scattering amplitude. 



APPENDIX E: APPLICABILITY OF THE TRANSPORT EQUATION 

The applicability of transport equations is rooted in Boltzmann's approximation of "molecular chaos" . A rigorous 
formulation of quantum field-theory (multiparticle quantum mechanics) in a context appropriate to study statistical 
mechanics can well start with the Wigner function, and deduce from it the Bogoliubov-Born-Grcen-Kirkwood-Yvon 
hierarchy of equations for multiparticle distribution functions. Decoupling the lowest order equation for the one- 
particle distribution function /(x, p, t) from the rest is an approximation of low density; two successive collisions of 
the same quantum must be uncorrelated. 

This condition is tantamount to stating that the mean free path (dependent on density and average cross section) 



A oc 



n(T)a 

is much smaller than the reach of the interaction, that can be expressed as the scattering length at low energies, or 
more generally as the square root of the averaged cross-section 



1 



>> Va 



n(T)a 
n{T)a 3/2 « 1 



(El) 



To satisfy this relation the interaction does not necessarily need to be very weak. As long as the free Bose density 
is a reasonable approximation (the system remains in a gaslike phase, as is the case for the LcrM), one can accept 
moderately strong interactions. Simple criteria for the strength of the interaction are to examine the ratio of the leading 
order (LO) and next-to-leading order (NLO) derivative expansion of the squared scattering amplitude Eq. (|B8[) 



LO 



NLO I 



2s 



2sN 2 



T7T ~ oo-a..a lo S -75 



32n 2 v 2 



(E2) 
(E3) 



or the ratio of the LO squared amplitude to the total square amplitude. The two inverse ratios are plotted in Fig. [T3] 
They clearly demonstrate that the interactions are indeed strong, thus requiring our large N resummation. We have 
however seen in Sec. IHII that the entropy density (and hence very likely the particle density) is close to the value in a 
free Bose gas. In Fig. [14] wc then plot the maximum possible value of the total cross-section, reached at the peak of 
the sigma resonance 



cm(s) 



32ir{v 2 ) 2 r 



1 - 



Mi 



' sN 
32-kv' 2 



(E4) 
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FIG. 14: Left panel: Typical cross section in the Linear Sigma model. We have taken an approximate average for Mandelstam 
s given by s ~ 10T 2 at a typical temperature T (with m T (0) = 0). Right panel: Particle number density for a free Bose gas in 
the same temperature range. 
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FIG. 15: The product a 3//2 n(T), significantly smaller than 1 as seen in this graph, shows that the Boltzmann-Uehling-Uhlenbeck 
equation is applicable in the range of temperatures that we have studied. 



and this particle number density 



n(T) = I Andpp 2 _ (E5) 

„/Vp 2 +™x( t ) _ 1 



Finally we construct the product n(T)(T 3//2 and plot it in Fig. 1151 The value of this parameter is seen to be of 
order 0.1 in our entire temperature range. Its actual value will be even smaller with an averaged cross-section (we 
are using its maximum possible value at a given temperature). Therefore there is no reason to question the validity 
of Boltzmann's equation in the model. 



APPENDIX F: SATURATION OF THE POLE MASS 



To further understand why the minimum of rj/s occurs just before the phase transition, we now perform a scaling 
analysis. For this, let us note that dimensionally and as a very rough average (not comparable with the quality 
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FIG. 16: The trailing tail of the a resonance occurs just before the phase transition temperature T c . Because the cross-section 
drops rapidly, the viscosity raises rapidly (faster than the entropy density s k T 3 ) and rj/s grows slightly before the phase 
transition. 



calculations presented in the main text) the scaling laws are 

S ~T 3 . (Fl) 

a 

From the cross section in Fig. [16] one can see that the maximum of the ittt cross section (sigma pole) occurs somewhat 
before the phase transition. The trailing tail of the resonance is well fit (see figure) by a power-law 

* = A 2+k 7 L, (F2) 

where A is a constant with [mass dimension] = 1. Then, just before the minimum: 

r)/s ~ T k - 2 . (F3) 

The fit gives the value k = 2.527 ± 0.020 that means that, for temperatures just above its minimum, rj/s ~ T 527 ; 
that is, the KSS coefficient begins to grow somewhat slowly toward the phase transition. The same conclusion can be 
reached by analytically examining the scattering amplitude in Eq. (|B8[) . 

Next we examine the connection of this turning of the cross section with the order parameter, the sigma condensate. 
In Fig. [T7] we plot the cross section that we have just commented on and the cross section where we have artificially 




FIG. 17: 7T7T cross section at finite temperature in the L<tM with (right) and without (left) the variation of the order parameter 
(<r)(T) that affects the amplitude in Eq. (JB8J) through F(T). 
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FIG. 18: 77/s and its dependence on the parameter Mrq, the tree-level (not the pole mass) a mass at T=0. The minimum of 
77/s increases with Mrq, but does not exceed T c . 
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FIG. 19: Saturation of the temperature at which the cross section is maximum and the temperature at which 77/s is minimum. 
Both are related because 77/s oc 1/a as shown in both Boltzmann and Green-Kubo formalisms. The critical temperature is 
always above these temperatures for any value of Mho- 



removed the variation of the order parameter (<t(T)), and fixed it as a constant. As can be seen, the variation of 
the order parameter with T makes the resonance narrower (and only slightly shifted toward smaller energies). This 
plotted behavior can be followed analytically from Eq. (|B8[) . Therefore, the slight rise of 77/s with temperature before 
the phase transition is indeed a direct consequence of the decrease of the order parameter. 

The pole mass position is highly dependent on Mrq, a parameter controlling the propagator at tree level. Both the 
maximum of the cross section and the minimum of 77/s move with Mrq. Just as in Fig. [5] we plot the dependence of 
77/s with Mm in Fig. QU 

In this figure we show that the exact temperature of the minimum grows with Mro- From this fact, one could ask 
if there exists a value of Mrq for which the minimum of 77/s is located at T c or above it. However our analysis stands 
as presented; increasing the value of Mro we observe that the pole mass in the squared amplitude saturates. The 
minimum of 77/s saturates as well when Mho 00. In Figure [15] we plot the temperatures at which the maximum 
of the cross section is located, and the one at which the minimum of 77/s is reached. Both values tend to a constant 
always below T c when Mrq is arbitrarily increased. We conclude that the minimum of 77/s takes place always below 
T c and never at or above it. 
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